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and Gaussianity 



O 

o 

(N 



O 

o 

I 

o 



M. Sadcgh Movahcd^-^^^, F. Ghasemi^, Sohrab Rahvar^-^, M. Reza Rahimi Tabar^'^ 
^Department of Physics, Sharif University of Technology, P.O.Box 11365-9161, Tehran, Iran 
^ Institute for Studies in theoretical Physics and Mathematics, P.O.Box 19395-5531, Tehran, Iran 
^ Iran Space Agency, P.O.Box 199799-43 13, Tehran, Iran 
*CNRS UMR 6202, Observatoire de la Cote d'Azur, BP 4229, 06304 Nice Cedex 4, France 

We investigate the statistical isotropy and Gaussianity of temperature fluctuations of Cosmic Mi- 
crowave Background radiation (CMB) data from Wilkinson Microwave Anisotropy Probe survey, 
using the Muhifractal Detrended Fluctuation Analysis, Rescaled Range and Scaled Windowed Vari- 
ance methods. These methods verify that there is no evidence for violation of statistical isotropy 
in CMB data. The multifractal detrended fluctuation analysis shows that Cosmic Microwave Back- 
ground fluctuations has a long range correlation function with a multifractal behavior. By comparing 
our analysis with the artificial shuffled and surrogate series of CMB data, we conclude that the mul- 
tifractahty nature of temperature fluctuation of CMB is mainly due to the long-range correlations 
and the map is consistent with a Gaussian distribution. 
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I. INTRODUCTION 

The Wilkinson Microwave Anisotropy Probe ( WMAP) 
mission is designed to determine the geometry, matter 
content, and the evolution of the universe. It is shown 
that the universe is geometrically flat and dark energy 
at the present time has the dominant contribution in 
the matter content of universe, causes the universe to 
accelerate rather than decelerate [1-4]. The statistical 
properties of the Cosmic Microwave Background radia- 
tion (CMB) data can be a unique tool to identify the the 
parameters of standard model of cosmology [5] . The re- 
sults of WMAP are a milestone in CMB anisotropy mea- 
surements since it combines high angular resolution, high 
sensitivity, with 'full' sky coverage allowed by a space 
mission. The frequency coverage of WMAP experiment 
allows CMB sky maps to be foreground cleaned up to 
I - 100 [4]. 

After the first year of WMAP data, the statistical 
isotropy and Gaussianity of the CMB anisotropy have 
attracted considerable attention. Tantalizing evidence of 
statistical isotropy breakdown and non-Gaussianity have 
mounted in the WMAP first year sky maps, using a va- 
riety of different methods. Since the observed CMB sky 
is a single realization of the underlying correlation, the 
detection of statistical isotropy violation or correlation 
patterns pose a great observational challenge. In order 
to extract more information from the rich source of in- 
formation provided by present (and future) CMB maps, 
it is important to design as many independent statisti- 
cal methods as possible to study deviations from stan- 
dard statistics, such as statistical isotropy. Since statis- 
tical isotropy can be violated in many different ways. 



various statistical methods can come to different con- 
clusions. Each method by design is more sensitive to 
a special kind of statistical isotropy violation. Recently 
bipolar power spectrum has shown no strong evidence 
of statistical isotropy violation [6,7], but analysis of the 
distribution of extrema in WMAP sky maps has indi- 
cated non-Gaussianity and to some extent, violation of 
statistical isotropy [8] . The Gaussianity of primordial 
fluctuations is a key assumption of modern cosmology, 
motivated by simple models of inflation [9-11]. The sta- 
tistical properties of the primordial fluctuations gener- 
ated by the inflationary cosmology arc closely related to 
the cosmic microwave background radiation anisotropy 
and a measurement of non-Gaussianity is a direct test 
of the inflation paradigm. If CMB anisotropy is Gaus- 
sian, then the angular power spectrum fully specifies its 
statistical properties. Many authors have searched for 
non-Gaussian signatures in CMB data using peak dis- 
tributions [12,13], the genus curve [14,15], peak corre- 
lations [16] and global Minkowski functionals methods 
[17]. More recently, the techniques used for the detec- 
tion of non-Gaussianity in hydrodynamic turbulence was 
applied for CMB data [18,19]. Moreover the Gaussian- 
ity of CMB in different angular scales have been tested 
[16,20-43]. Most of the previous works only tested the 
consistency between the CMB data and simulated Gaus- 
sian realizations and so far they found no significant ev- 
idence for cosmological non-Gaussianity. 

In this paper we would like to characterize the complex 
behavior of CMB through the computation of the fluctu- 
ation parameters - scaling exponents - which quantifies 
the correlation exponents and multifractality of the data. 
We use certain fractal analysis approaches such as, the 



Multifractal Detrended Fluctuation Analysis (MF-DFA), 
Rescaled Range Analysis (R/S) and Scaled Windowed 
Variance (SWV) to analysis the data set. Using the new 
approaches, we will test the statistical isotropy and Gaus- 
sianity of temperature fluctuations at the last scattering 
surface. The MF-DFA method shows that Cosmic Mi- 
crowave Background fluctuations has a long range cor- 
relation function with multifractal behavior. Comparing 
the MF-DFA results of the original temperature fluctua- 
tions to those for shuffled and surrogate series, we con- 
clude that the multifractality nature of CMB is mainly 
due to long-range correlations and the map is consistent 
with Gaussian distribution. Applied methods of MF- 
DFA, R/S and SWV, show that WMAP data is a sta- 
tistically isotropic data set or fractional Gaussian noise. 
The value of Hurst exponent guarantees that there is no 
evidence for violation of statistical isotropy in the CMB 
anisotropy map considered here. 

This paper is organized as follows: In Section II we 
briefly describe the MF-DFA method and show that the 
scaling exponents determined via the MF-DFA method 
are identical to those obtained by the standard multifrac- 
tal formalism based on partition functions. Also briefly 
introduce the Rescaled Range Analysis and Scaled Win- 
dowed Variance method. The fractal analysis of the tem- 
perature fluctuations are presented in Section III. In Sec- 
tion IV, we compare the multifractal behavior of the orig- 
inal data with that of shuffled and surrogate series and 
show that the multifractality is mainly due to the long- 
range correlations. In Section V we investigate the sta- 
tistical isotropy and Gaussianity of temperature fluctua- 
tions using the MF-DFA, R/S and SWV results. Section 
VI closes with a discussion of the present results. 



II. FRACTAL ANALYSIS METHODS 

In this section we introduce three methods to investi- 
gate the fractal properties of stochastic processes. 

A. Multifractal Detrended Fluctuation Analysis 

The MF-DFA methods are the modifled version of de- 
trended fluctuation analysis (DFA) to detect multifractal 
properties of series. The detrended fluctuation analysis 
(DFA) method introduced by Peng et al. [44] has became 
a widely-used technique for determination of (mono)- 
fractal scaling properties and the detection of long-range 
correlations in noisy non-stationary time series [44-48]. 
It has successfully been applied to diverse fields such 
as DNA sequences [44,49], heart rate dynamics [50,51], 
neuron spiking [52], human gait [53], long-time weather 
records [54], cloud structure [55], geology [56], ethnology 



[57], economical time series [58], solid state physics [59] 
and solar physics [60]. 

The simplest type of the multifractal analysis is based 
upon the standard partition function multifractal formal- 
ism, which has been developed for the multifractal char- 
acterization of normalized, isotropic (stationary) mea- 
surements [61-64]. Unfortunately, this standard formal- 
ism does not give correct results for non-isotropic angular 
(non-stationary time) series that are affected by trends or 
that cannot be normalized. Thus, in the early 1990s an 
improved multifractal formalism has been developed, the 
wavelet transform modulus maxima (WTMM) method 
[65] , which is based on the wavelet analysis and involves 
tracing the maxima lines in the continuous wavelet trans- 
form over all scales. The other method, the multifractal 
detrended fluctuation analysis (MF-DFA), is based on 
the identification of scaling of the gth-order moments de- 
pending on the signal length and is a generalization of 
the standard DFA using only the second moment q = 2. 

The MF-DFA does not require the modulus maxima 
procedure in contrast to the WTMM method, and hence 
does not require more effort in programming and comput- 
ing than the conventional DFA. On the other hand, of- 
ten experimental data are affected by non-isotropic (non- 
stationarities) like trends, which have to be well distin- 
guished from the intrinsic fluctuations of the system in 
order to find the correct scaling behavior of the fluctua- 
tions. In addition very often we do not know the reasons 
for underlying trends in collected data and even worse 
we do not know the scales of the underlying trends, also, 
usually the available recorded data is small. For the reli- 
able detection of correlations, it is essential to distinguish 
trends from the intrinsic fluctuations in the data. Non- 
detrending methods work well if the records are long and 
do not involve trends. But if trends are present in the 
data, they might give wrong results. Detrended fluctua- 
tion analysis (DFA) is a well-established method for de- 
termining the scaling behavior of noisy data in the pres- 
ence of trends without knowing their origin and shape 
[44,50,66-68] 

The modified multifractal DFA (MF-DFA) procedure 
consists of five steps. The first three steps are essentially 
identical to the conventional DFA procedure (see e. g. 
[44-48]). In our case which is studying the temperature 
fiuctuations of CMB, we take the temperature series with 
the size of N and follow the steps as follows: 
• Step 1: Determine the "profile" 

y(7.) ^ ^ [r(nO - (T)] , s = l,...,iV, (1) 

1=1 

where Tifii) is the temperature of CMB map, lii is the 
unique vector pointing CMB and 7s = arccos(ni.ns) is 
the size of segment of CMB that we are calculating the 
series. Subtraction of the mean (T) is not compulsory, 
since it would be eliminated by the later detrending in 



the third step. 

• Step 2: Divide the profile Y{-fs) into N^^ = int(A^/s) 
non-overlapping segments of equal angular lengths 7^. 

• Step 3: Calculate the local trend for each of the N^^ 
segments by fitting a polynomial function to y(7s). The 
variance between l^(7s) and the function of the best fit 
for each segment v ^ 1, • ■ • , -/V^^ , is as follows: 

F2(7„ v)^-Y^ {Y\{v 1)7. + 7.] - vMf , (2) 

i=\ 

A Linear, quadratic, cubic, or higher order polynomials 
can be used in the fitting procedure (conventionally called 
DFAl, DFA2, DFA3, . . .) [44,46,47,51]. 

• Step 4- Average over all segments to obtain the g-th 
order fluctuation function, defined as: 

FM^l^^Y.[FHjs,,^)f'\ , (3) 

where, in general, the index variable q can take any real 
value except zero. For q = 2, the standard DFA pro- 
cedure is retrieved. Generally we are interested in how 
the generalized q dependent fluctuation functions Fq{'^s) 
depend on the angular scale 7^ for different values of q. 
Hence, we must repeat steps 2, 3 and 4 for several an- 
gular scales 7s. It is apparent that ^"^(7^) will increase 
with the increasing of 7s . 

• Step 5: Determine the scaling behavior of the fluctua- 
tion functions by analyzing log- log plots of ^^5(7^) versus 
7s for each value of q. If the series Tirii) are long-range 
power-law correlated, Fq{"fs) increases, for large values of 
7s, as a power-law 

(4) 

In general, the exponent h{q) may depend on q. For 
isotropic fluctuations, < h{2) < 1.0 and h{2) is iden- 
tical to the well-known Hurst exponent {h{2) ~ H) 
[44,45,61]. In the absence of statistical isotropy the corre- 
sponding scaling exponent of i^g(7s) is larger than unity 
h{2) > 1.0 and its relation to the Hurst exponent is 
H = h{q = 2) - \ [44,60,69]. Thus, one can caU the 
function h{q) as the generalized Hurst exponent. 

For monofractal fluctuations, h{q) is independent of g, 
since the scaling behavior of the variances (7s , v) is 
identical for all segments of and the averaging proce- 
dure in Eq. (3) will just give this identical scahng be- 
havior for all values of g. If we consider positive values 
of g, the segments u with large variance F'^{'^s,v) (i. e. 
large deviations from the corresponding fit) will dominate 
the average Fq{"fs)- Thus, for positive values of q, h{q) 
describes the scaling behavior of the segments with large 
fluctuations. Usually the large fluctuations are character- 
ized by a smaller sealing exponent h{(i) for multifractal 
series [70]. On the contrary, for negative values of g, the 



segments v with small variance F'^ (7s , v) will dominate 
the average -F'q(7s). Hence, for negative values of q^ h{q) 
describes the scaling behavior of the segments with small 
fluctuations [70]. 

1. Relation to standard multifractal analysis 

For a isotropic series the multifractal scaling expo- 
nents h{q) defined in Eq. (4) are directly related to the 
scaling exponents T{q) defined by the standard partition 
function-based multifractal formalism as shown below. 
Suppose that the series T{hi) of length N is an isotropic 
sequence. Then the detrending procedure in step 3 of 
the MF-DFA method is not required, since no trend has 
to be eliminated. Thus, the DFA can be replaced by the 
standard Fluctuation Analysis (FA) with the definition 
of variance for each segment ly, ly = I, . . . , N^^ as follows: 

^2^(7., 1^) ^ [Y{,yjs) - Y{{u - l)7s)]^ (5) 

Inserting this simplified definition into Eq. (3) and using 
Eq. (4), we obtain 

-^Y.\Y{v^,)-Y{{v-~l)^,)A ^7s'^"- (6) 

u—l J 

In order to relate also to the standard textbook box 
counting formalism [61-64], we employ the definition 
of the profile in Eq. (1). It is evident that the term 
y(i/7s) — Y{{v — l)7s) in Eq. (6) is identical to the sum 
of the numbers Tihi) within each segment v of size 7s. 
This sum is known as the box probability p^^ (y) in the 
standard multifractal formalism for T(hi) 

us 

p,M^ E T{n,)^Y{,y^,)-Y{{,y~l)js), (7) 

The sealing exponent T(g) is usually defined via the 
partition function Zq{'^s) 

Z,(7s)^ElP'r,Mr~7j''\ (8) 

where q is a real parameter as in the MF-DFA method, 
discussed above. Using Eq. (7), we see that Eq. (8) is 
identical to Eq. (6), and obtain analytically the relation 
between the two sets of multifractal scaling exponents 

T{q) = qh{q) - 1. (9) 

Thus, we observe that h{q) defined in Eq. (4) for the 
MF-DFA is directly related to the classical multifractal 
scaling exponents T{q). Note that ti{q) is different from 
the generalized multifractal dimensions 



(10) 



that are used instead of t((j) in some papers. While h{q) 
is independent of q for a monofractal time series, D{q) 
depends on q in this case. Another way to characterize a 
multifractal series is the singularity spectrum f{a), that 
is related to r(g) via a Legendre transform [61,63] 



T'{q) and f{a)=qa-T{q). 



(11) 



here, a is the singularity strength or Holder exponent, 
while f{a) denotes the dimension of the subset of the 
series that is characterized by a. Using Eq. (9), we can 
directly relate a and f{a) to h{q) 

a ^ h{q) + qh'{q) and f{a) = q[a - h{q)] + 1. (12) 

A Holder exponent denotes mono-fractality, while in the 
multifractal case, the different parts of the structure are 
characterized by different values of a, leading to the ex- 
istence of the spectrum /(a). 



(15) 



where r and H arc time span and Hurst exponent, re- 
spectively. For temperature fluctuations on the last scat- 
tering surface, R and S are defined as 



i?(7,) = Max{y(7,)} - Min{y(7,)}, 

1/2 



(16) 



^(^«) ^{hi ~ ^^^]') ' ^ = 1' • ■ • ' (17) 

where F(7s) is defined according to Eq. (1). The scal- 
ing behavior of the fluctuation function is determined by 
analyzing log-log plot oi R/ S versus 7s as 



R/S ^ 7f ■ 



(18) 



III. FRACTAL ANALYSIS OF COSMIC 
MICROWAVE BACKGROUND RADIATION 
DATA 



B. Scaled Windowed Variance Analysis 

The Scaled Windowed Variance analysis was developed 
by Cannon et al. (1997) [69]. The profile of temperature, 
y(7s), is divided into intervals of angular length scale 7s. 
Then the standard deviation is calculated within each 
interval using the following relation 

/ s 

SWV(7s) = f - E[^(^') - '^Y{ls))f\ ■ (13) 

The average standard deviation of all angular intervals 
of length 7s is computed. This computation is repeated 
over all possible interval lengths. The scaled windowed 
variance is related to 7s by a power law 

SWV - 7s^. (14) 



C. Rescaled Range Analysis 

Hurst developed the Rescaled Range analysis, a statis- 
tical method to analyze long records of natural phenom- 
ena [69,71]. There are two factors used in this analysis: 
firstly the range i?, this is the difference between the min- 
imum and maximum 'accumulated' values or cumulative 
sum of X{t,T) of the natural phenomenon at discrete 
time t over a time span r, and secondly the standard 
deviation 5, estimated from the observed values X{t,T). 
Hurst found that the ratio R/S is very well described for 
a large number of natural phenomena by the following 
empirical relation 



As mentioned in section II, a spurious of correlations 
may be detected in the absence of statistical isotropy ( 
if series is non-stationarity), so direct calculation of cor- 
relation functions, fractal dimensions etc., may not give 
reliable results [44,45,61,69,70]. The simplest way to ver- 
ify the statistical isotropy of the temperature fluctuations 
on the last scattering surface is by measuring the stabil- 
ity of the variance of the temperature in various sizes of 
windows. Figure 1 shows the standard deviation of the 
temperature verses the angular size 7s of the window. 
Here we have a saturation for the standard derivation 
in the large angular scale which reflects the long-range 
correlation behavior of the temperature fluctuations of 
CMB [44,45,61,69,70]. 

We use the MF-DFAl method for analyzing the tem- 
perature fluctuations on CMB. Following the procedure 
for MF-DFA analysis as described in the last section we 
obtain -F'q(7s) as a function of angular scale 7s. For each 
index of g the generalized Hurst exponents h{q) in Eq. (4) 
can be found by analyzing log-log plots of Fq^-jg) versus 
7s. Figure 2 shows the Hurst exponent in terms of q for 
MF-DFAl analysis. Variation of Hurst exponent in terms 
of q shows the multifractal behavior of temperature fluc- 
tuation of CMB. More over we have different values of 
the slope of classical multifractal scaling exponent ,t((7), 
for g < and q > 0, which indicates that CMB has a 
multifractal structure. 

For positive and negative values of q, we obtain T{q) 
with the slopes of 0.79 ± 0.03 and 1.00 ± 0.03, respec- 
tively. According to the relation between the Hurst expo- 
nent and MF-DFA exponent, we obtain H = 0.94 ± 0.01 
which means that temperature fluctuation is an isotropic 
process with long-range correlation [70] . The variation of 




FIG. 1. Behavior of standard deviation of temperature 
fluctuations of CMB as a function of angular scale, 7^. 



singularity spectrum f{a), Eq. (11) also is shown in Fig- 
ure 2. The values of derived quantities from MF-DFAl 
method, are given in Table I. 

We also use the Scaled Windowed Variance (SWV) 
analysis [69,71] to determine the Hurst exponent for 
CMB via Eq. (14). Figure 3 shows log-log plot of Scaled 
Windowed Variance of CMB fluctuations as a function 
of angular scale, 7^ which results H = 0.95 ± 0.02. Fi- 
nally we apply the Rescalcd Range analysis (R/S), to 
determine the Hurst exponent of CMB fluctuations. Ac- 
cording to Figure 4, the value of Hurst exponent obtain 
as 0.95 ± 0.02, which is in agreement with the two previ- 
ous results. 
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FIG. 2. The q dependence of the generalized Hurst expo- 
nent h{q), the corresponding r(q) and singularity spectrum 
/(a) are shown in the upper to lower panel respectively for 
temperature fluctuation series. 



IV. COMPARISON OF MULTIFRACTAL 
NATURE FOR ORIGINAL, SHUFFLED AND 
SURROGATE CMB DATA 

In this section we are interested in to determine the 
source of multifractality and test the Gaussianity of CMB 
data. In general, two different types of multifractality in 
certain data set can be distinguished: (i) Multifractality 
due to a fatness of probability density function (PDF) 
and (ii) Multifractality due to different correlations in 
small and large scale fluctuations. In the first case the 
multifractality cannot be removed by shuffling the data 
set, while in the second case data may have a PDF with 
finite moments, e. g. a Gaussian distribution and the cor- 
responding shuffled series will exhibit mono-fractal scal- 
ing, since all long-range correlations are destroyed by the 
shuffling procedure. 

If we have both kinds of multifractality in data, the 
shuffled series will show weaker multifractality than the 
original series. The easiest way to distinguish the type of 
multifractality, is by analyzing the corresponding shuffled 
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FIG. 3. Behavior of Scaled Windowed Variance (SWV) of 
CMB fluctuations as a function of angular scale, 73 in log-log 
plot. 
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FIG. 4. The log-log plot of Rescaled Range (R/S) versus 
angular scale, 7^ for temperature fluctuation at last scattering 
surface. 

TABLE I. The values of the Hurst, multifractal scaling and 
generalized multifractal exponents for q — 2.0, for original, 
surrogate and shuffled of temperature fluctuation series ob- 
tained by MF-DFAI. 



Data 


H 


r 


D 


CMB 


0.94 ±0.01 


0.88 ±0.02 


0.88 ±0.02 


Surrogate 


0.88 ±0.01 


0.76 ±0.02 


0.76 ±0.02 


Shufficd 


0.50 ±0.001 


0.002 ± 0.002 


0.002 ± 0.002 



and surrogate series. The shuffling of data set destroys 
the long range correlation, therefore if the niultifractal- 
ity only belongs to the long-range correlation, we should 
find /ishufl?) = 0.5 [70]. The multifractality nature due 
to the fatness of the PDF signals is not affected by the 
shuffling procedure. On the other hand, to determine the 
multifractality due to the broadness of PDF by the sur- 
rogating method, the phase of discrete fourier transform 
(DFT) coefficients of CMB data are replaced with a set 
of pseudo independent uniform distribution within the 
range of (— tt, tt). The correlations in the surrogate series 
do not change, while the probability function changes to 
the Gaussian distribution. If multifractality in the data 
set is due to a broadness of PDF, h{q) obtained by the 
surrogate method will be independent of q. If both kinds 
of multifractality are present in CMB data, the shuffled 
and surrogate series will show weaker multifractality than 
the original one. 

To check the nature of multifractality, we compare the 
fluctuation fimction Fq{'^s), for the CMB map with the 
result of the corresponding shuffled, F|^"^(7s) and sur- 
rogated data F™''(7s). Differences between these two 
fluctuation functions with the original one can indicate 



the presence of long range correlations or broadness of 
probability density function in the original data. These 
differences can be obtained by the ratio F^(7s)/-F|^"'(7s) 
and Fq{'-fs)/F™'^{'ys) as a function of 7s [70]. Since 
the anomalous scaling due to a broad probability den- 
sity affects Fq{'fs) and -F|'^"^(7s) in the same way, only 
multifractality due to correlations will be observed in 
Pqils)/ Fq^^^ {'js)- The scaling behavior of these ratios 
are: 

^^,(7.)/i^f "'(7s) ^ 7^^^-"=""^^'^^ = l^""^'''- (19) 

Fqijs)/F^^^r.) ^ 7M?)-'^s„(<,) ^ ^/.Pop(,). (20) 

If only fatness of the PDF be responsible for the multi- 
fractality, one should find h{q) = /ishuf('?) and ft.corr('z) = 
0. On the other hand, deviations from ft.corr('?) = in- 
dicates the presence of correlations, and q dependence of 
hcoi-riq) indicates that multifractality is due to the long- 
rage correlation. 

If the multifractal behavior of CMB is made by 
the broadness of PDF and long-range correlation, both 
/ishuf('?) and hsui {q) will depend on q. The q dependence 
of the generalized exponent h(q) for original, surrogate 
and shuffled CMB data are shown in Figures 5 which 
shows that the multifractality nature of temperature fluc- 
tuation is almost due to the long — range correlation. 
The absolute value of hconiq) is larger than /ipdf(<7) 
which indicates that the multifractality due to the fat- 
ness is much weaker than the multifractality due to the 
correlation. The deviation of /isur(9) and hshuiil) from 
h(q) can be determined by using as follows: 

^^a{q,r + a4q,r' ^''^ 

where the symbol " o " represents for "surroated" 
and "shuffled" series. The value of reduced chi-square 
xio = xl/-^ is the number of degree of freedom) for 
shuffled and surrogate time series are 21599.11, 313.43, 
respectively. 

Singularity spectrum f{a) of the surrogate series is 
almost similar to the original temperature fluctuations, 
while in the case of shuffled series we have a narrower sin- 
gularity spectrum of Aashuf = "(g^in) - a{qmax) - 0.02 
compare to that of surrogate Ackgur — 0.09 and original 
data Aa ~ 0.12. The small value of Aaghuf compare to 
the two other series shows that the multifractality in the 
shuffled CMB map has almost been destroyed [72]. Fig- 
ures 6 and 7 show the MF-DFAl results for surrogate and 
shuffled temperature fluctuation series. Comparing the 
MF-DFA results of the data set to those for shuffled and 
surrogate series, we conclude that the multifractal na- 
ture of temperature fluctuations in CMB data is almost 
due to the long — range correlations and the probabil- 
ity distribution function of CMB is almost Gaussian (see 
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FIG. 5. Generalized Hurst exponent as a function of q for 
original, surrogate and shuffled CMB data. 



also the section V). The values of the Hurst exponent 
if, multifractal scaUng T{q = 2), generalized multifractal 
dimension D{q = 2) of the original, shuffled and surro- 
gate data of CMB obtained with MF-DFAl method are 
reported in Table I. 



V. STATISTICAL ISOTROPY AND 
GAUSSIANITY OF CMB ANISOTROPY 
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FIG. 6. The q dependence of the generalized Hurst expo- 
nent the corresponding r((jr) and singularity spectrum 
/(q) are shown in the upper to lower panel respectively for 
surrogate temperature fluctuation series. 



In this section we examine the Gaussianity and Statis- 
tical isotropy of temperature fluctuations in CMB data 
tlirough MF-DFA, RjS and SWV methods. 

Standard models of inflationary cosmology predict a 
Gaussian distribution of CMB temperature fluctuations. 
However, along witli standard inflationary models, there 
exist theories such as inflationary scenarios with two or 
more scalar fields which predict a non-Gaussian primor- 
dial fluctuations [73-76] . Another possibility of deviation 
from the Gaussianity is the manipulation of the CMB 
signals after the recombination due to subsequent weak 
gravitational lensing [77,78] and various foregrounds like 
dust emission, synchrotron radiation, or unresolved point 
sources [79] . One should also take into account the addi- 
tional instrumental noise in the observational data [80]. 
One of the advantages of Gaussian field is that the two 
point correlation function C(n, n') = {T{n)T{n')) will be 
sufficient for fully specify the statistical properties of the 
field. In the case of non-Gaussianity, one must take the 
higher moments of the field into account in order to fully 
describe the whole statistics. So, studying the Gaussian 
nature of the signals or detecting some distinctive non- 
Gaussianity is an important issue in CMB physics. 

As we have seen in Sec. IV, the MF-DFAl results of the 
shuffled series shows almost constant h{q) ~ 0.5 which 
indicates that h{q) is independent of the q or in another 
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FIG. 7. The q dependence of the generalized Hurst expo- 
nent h{q), the corresponding r(g) and singularity spectrum 
f{a) are shown in the upper to lower panel respectively for 
shuffled temperature fluctuation series. 



word the shuffled series is a monofractal structure. This 
means that tlic multifractality of the original data is due 
to the long-range correlations not due to a broadness of 
PDF of temperature fluctuations (Figure 5). Further- 
more small differences in the generalized Hurst exponent 
between the original and surrogated data (Figure 5) and 
similarity of the singularity spectrum, /(a), in those two 
series, show that the probability density function of tem- 
perature fluctuations has well confidence with the Gaus- 
sian distribution which is in agreement with the recent 
results [6,7,18,19] 

The second important question about the CMB data 
is its statistical isotropy. One of the main motivations for 
anisotropy of the universe comes from the global topology 
of the universe. The idea of a multi-connected topology of 
the universe comes back to Schwarzschild's work [81] be- 
fore the Einstein's first seminal paper on cosmology [82] . 
If the universe is assumed to be simply-connected, the 
cosmological principle implies that it should be globally 
isotropic to any observer. In the case of multi-connected 
universe while the Einstein gravity still hold in this space 
but the global structure of universe at large scales will 
be more complicated. One of the most promising sig- 
natures may be on the CMB as the existence of large 
scale anisotropies with a repeating patterns. In addition 
to cosmological mechanisms, instrumental and environ- 
mental effects in observations of CMB can also easily 
produce spurious breakdown of statistical isotropy. Here 
we examine the anisotropy of CMB with the MF-DFA 
analysis. 

The statistical isotropy means that the statistical prop- 
erties of CMB (e.g. n-point correlations) to be invariant 
under the Eulerian rotation. In another word: 

C(ni,n2) = C(ni-n2)=C(7), (22) 

where 7 = arccos(rii •^2) is the angle between hi and ^2- 
It is then convenient to expand it in terms of Legendre 
polynomials, 

-| 00 

^(^) = ^ E(2^ + mmcosi), (23) 
1=2 

where Ci is the widely used angular power spectrum. The 
summation over I starts from 2 because the / = term is 
the monopole which in the case of statistical isotropy the 
monopole is constant, and it can be subtracted out. The 
dipole / = 1 is due to the local motion of the observer 
with respect to the last scattering surface and can be 
subtracted out as well. In order to extract the statistical 
property of CMB we need an ensemble maps from the 
CMB to calculate the correlation functions. In reality 
we have only one map from the CMB, but if we assume 
the statistical isotropy we can calculate the correlation 
functions moving on overall data. 

As we mentioned in section II, for a process with sta- 
tionarity in time or isotropic process, the Hurst exponent 



is always less than unity [44,69]. According to the MF- 
DAFl, SWV and R/S methods we obtained the value 
of Hurst exponent for WMAP data as = 0.94 ± 0.02. 
These analysis with the Hurst exponent {Q < H < 1.0) 
show that there is no evidence for violation of statistical 
isotropy of WMAP data. 

VI. CONCLUSION 

In this work we used three independent approach to in- 
vestigate the statistical properties of WMAP data. The 
applied methods of Multifractal Detrended Fluctuation 
Analysis, Rescaled Range Analysis and Scaled Windowed 
Variance methods show that the WMAP data is a statis- 
tically isotropic data set or a fractional Gaussian noise. 
The value of Hurst exponent, H ~ 0.94, guarantees that 
there is no evidence for violation of statistical isotropy 
in the CMB map. The MF-DFA method allows us to 
determine the multifractal characterization of this map. 
The q dependence of generalized Hurst h{q) and classical 
multifractal scaling r(g) exponents, show that the tem- 
perature fluctuations on the last scattering surface has 
multifractal behavior which is agreement to the extended 
self-similarity method [18]. We obtained the generalized 
multifractal dimension of CMB data D ^ 0.88. Accord- 
ing to a small difference between generalized Hurst ex- 
ponent and singularity spectrum of the original data set 
with the surrogate one's, H — Hsm- = 0.06 ± 0.01 and 
Aa — Aofsur = 0.03, we found that there is no crucial 
evidence for non-Gaussianity of the CMB map. Com- 
parison of the generalized Hurst exponent of the original 
data with the shuffled and surrogate series showed that 
the nature of the multifractality of CMB is mainly due 
to the long-range correlation , rather than the broadness 
of the probability density function. 
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